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We cast the infinitesimal generator of the symmetric simple exclusion process on a graph as a quantum spin- 
2 ferromagnetic Heisenberg model and show that its eigenspectrum can be obtained by elementary techniques 
on the complete, complete bipartite, and related multipartite graphs. Some of the resulting infinitesimal gen- 
erators are formally identical to homogeneous as well as mixed higher spins models. The degeneracies of the 
eigenspectra are described in detail, and we give a neat derivation of the outer multiplicities appearing in the 
Clebsch-Gordan series for arbitrary spin-s representations of the SU(2) not easily found elsewhere. We mention 
in passing how our results fit within the related questions of a ferromagnetic ordering of energy levels and a 
conjecture according to which the spectral gaps of the random walk and the interchange process on finite simple 
graphs must be equal. 
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I. INTRODUCTION 

Exclusion processes, together with the contact process and 
the Glauber-Ising model, are one of the most fundamental 
models in the field of nonequilibrium interacting particle sys- 
tems [1]. In physics, exclusion processes are the simplest 
models that furnish nontrivial results on a number of basic 
issues, such as the relaxation dynamics of an interacting gas 
towards the thermodynamic equilibrium or the dynamics of 
shock waves in discrete models for inviscid fluids [2]. In the 
one-dimensional linear chain, simple exclusion processes, ei- 
ther symmetric or asymmetric, under periodic or more gen- 
eral open boundary conditions, have been analyzed and their 
relationship with other models of interest spanned a wealth of 
mathematical physics during the last two decades [3-9]. 

The investigation of exclusion processes on general graphs, 
however, has been timid in the physics literature, except per- 
haps for the study of some multi-lane traffic models that in- 
clude junctions and bifurcations [10], and most studies do not 
go beyond the level of simple mean field analyses or numeri- 
cal simulations. In the mathematical literature, otherwise, the 
study of random walks and exclusion processes on graphs is a 
relatively hot topic connected with deep results in probability, 
group theory, harmonic analysis, and combinatorics [11-13]. 
Unfortunatelly, this mathematical literature is inaccessible to 
most practicing physicists, even the well prepared ones, with 
possibly useful results buried beneath loads of advanced pre- 
requisites, hardcore formalism, and subtle rationale. 

The objective of this article is to show that the infinitesimal 
generator of the symmetric simple exclusion process (SSEP) 
on the complete, complete bipartite, and closely related graphs 
can be solved by elementary techniques that belong in the 
toolbox of every trained physicist. We believe that the ex- 
plicit calculations presented here simplify the understanding 
of the models and also open some interesting perspectives. 



II. THE SSEP ON A GRAPH 

Let G = (V,E) be a finite simple (without loops) undirected 
connected graph of order N with vertex set V = { 1 , . . . , N} and 
edge set E C V x V . To each i'eVwe attach a random variable 
<7[ taking values in {— 1,+1}. If a, = —1 we say that vertex 
i is empty and if Oi = +1 we say that vertex i is occupied by 
a particle. The state of the system is specified by the configu- 
ration a = (<j u ...,g n ) in Q. = {-l,+l} y . The SSEP(G) is 
the continuous -time Markov jump process that describes the 
transitions of a set of n itinerant particles, 1 ^ n ^ N, between 
the connected vertices of G. In the SSEP(G), each particle 
chooses, sequentially and at exponentially distributed times, 
one of its adjacent vertices to jump to provided the target ver- 
tex is empty, otherwise the jump attempt fails and the process 
continues. Clearly, when n = 1 we have the simple random 
walk on G. When n ^ 2, exclusion between particles comes 
into play and the process becomes more interesting. 

We introduce vector spaces in the description of the 
SSEP(G) by turning Q. into (C 2 ) 0V , <7 into |rj) = |oi) ® ••• 
® \ <Jn), and setting |0) = (J) and |1) = (i) to identify respec- 
tively an empty and an occupied vertex. A little reflection 
shows that within this vector space scenario the infinitesimal 
generator of the time evolution of the SSEP(G) can be written 
as (for detailed derivations see, e.g., [7-9]) 

# = £(i-%), (l) 

where i ~ j stands for pairs of connected vertices of G and fP, ; - 
is the operator that transposes the states of vertices and j, 

,<*{,-•• >°)V") = I"- >°)V ,<*,•••)• (2) 

As is well known, 2^, can be written in terms of Pauli spin 
matrices as 
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% } = 1(1 + o r oj) = 1(1 + of a] + a) a] + of a;). (3) 



2 



Inserting this fP, ; - in Eq. (1) gives 

#=*£(1-C7 r 3;). (4) 

We then see that is, to within a diagonal term, exactly the 
Hamiltonian of the isotropic Heisenberg spin-j quantum fer- 
romagnet over G [14]. The ground states of 9{ have eigen- 
value zero and correspond (under a probabilistic normaliza- 
tion) to the stationary states of the process. 

Operator (4) is positive semi-definite and the master equa- 
tion governing the time evolution of the probability den- 
sity P{o,t) of observing configuration a at instant t reads 
d,P(a,t) = —9-(P(<3,t). One is usually interested in the spec- 
tral gap of 9{, which is the inverse of the leading characteristic 
time scale of the process related with the time it takes to ap- 
proach the stationary state. Conservation of particles in the 
SSEP(G) implies that 9{ commutes with the total number of 
particles operator 

^=oB 1 + CJ r) = y+^ (5) 

where S z is the z-axis "polarization" operator. It follows 
that H is block-diagonal, H = ®„ Hi, with each block H 
acting on its respective invariant subspace Q.„ of dimension 
dim£2„ = f ) = N\/n\{N — n)\. The eigenspectrum of 9{ is 
also symmetric about n = N/2, because it commutes with the 
"spin flip" operator 

U = Y[of, (6) 

that transforms particles into holes and vice-versa, 
11\g\, ■ ■ ■ , <Jn) = \— Oir" > — Gn)> taking a state with n 
particles into a state with N — n particles. The eigenspectra in 
the sectors of n and N — n particles are thus identical. 

In what follows we investigate operator (4) on a couple of 
different graphs and show that some of them can be analyzed 
by elementary techniques. 

m. THE SSEP ON THE COMPLETE GRAPH 

In the complete graph Kn, every pair of distinct vertices is 
connected by a unique edge; see Fig. 1. For this graph, the 
infinitesimal generator (4) reads 

rt=\ E (l-S,-Bj). (7) 

Rearranging the summation in Eq. (7) and noticing that 
(<7, ) 2 = 3/ we eventually arrive at 

*=?(f+o-at<- <•> 

This 9i is but the Curie- Weiss version of the spin-j ferro- 
magnetic Heisenberg model without the overall multiplicative 



1 /N term usually included to keep the energy per spin an in- 
tensive quantity, since we are not doing any thermodynamics 
here [15]. In the basis simultaneously diagonal in the total 
spin squared operator 

/ 1 N \2 

with eigenvalues S(S + 1), S = 5 m ; n , S m { n + l,...,N/2, where 
Smin = or 1 /2 depending whether is even or odd, and in the 
total z-axis component S z defined in Eq. (5) with eigenvalues 
M = —S, —5 + 1 , . . . , +S, the eigenvalues of 9l read 

£ JV (5,M) = |(^ + l)-5(5+l), (10) 

i.e., En(S,M) = En(S), independent on M. The degeneracy 
of En(S) is given by gp/(S) = (25+ 1) x di/ 2 {N,S), where 
the factor 25+1 comes from the degeneracy in the S z values 
of rotationally invariant operators like !H, and the din(N,S) 
comes from the fact that there exists many possible combi- 
nations of the elementary spins summing up to a definite 
value of S. This last factor is given by the outer multiplicity of 
the irreducible representation (D^ appearing in the Clebsch- 
Gordan series 

N/2 

[2>(l/2)fW = d l/2 {N,S)rt s \ (11) 
and can be shown to be given by (cf. appendix) 

d ^ N ^={iN + S )-^N+S + l} (12) 

We see from Eqs. (10) and (12) that E N (S =N/2) = with 
a (N + l)-fold degeneracy. These values have a simple in- 
terpretation: the SSEP(G) has a zero eigenvalue on each of 
its jV + 1 sectors of total particle number n = N/2 + M = 
0, 1, . . . ,iV. That the stationary states of 9{ occur in the sec- 
tors of S — N/2 is just another statement of the well known 
fact that the ground states of ferromagnetic Heisenberg mod- 
els have maximum possible total S. The right eigenvectors 
corresponding to the zero eigenvalues are the stationary states 
of the process, explicitly given by 

KM) =(") ' I \i h ,i h) - (13) 

\ n / l<i'i</ 2 <- -<i n <N 

notice the probabilistic normalization of \<f$(n)), not the 
quantum-mechanical one. The summation in Eq. (13) runs 
over all combinations of the n particle positions z'i, i 2 , . . • , i n 
among the available vertices of the graph. 

For processes that conserve the total number of parti- 
cles like the SSEP(G), a basis diagonal in n is more use- 
ful. In the \S,M) basis, each invariant subspace Cl„ of 
fixed n = 0,1,..., N is spanned by the states with M = 
-N/2 + n fixed and \-N/2 + n\ ^ S < N/2, with the given 
\S,M = — N/2 + n) states within £2„ bearing their original 
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FIG. 1. Left: Complete graph A*g with n = 3 particles (black circles) 
occupying some of its vertices. Right: Complete bipartite graph A3 5 
occupied by n = 4 particles. 

multiplicity d\/2{N,S). This completely characterizes the 
eigenspectrum of 9-f in each of its invariant subspaces. 

Tables I and II illustrate the SSEP(X^) in the concrete case 
of N = 8. The eigenvalues of SSEPtXg) and their degenera- 
cies appear in Tab. I. The eigenspectrum in terms of the total 
number of particles appears in Tab. II and is clearly symmetric 
about n = N/2, as we anticipated in Sec. II. 

The spectral gap of H is given by the smallest nonzero 
eigenvalue of 9{, and is related with the characteristic time Taj 
it takes for the process to exponentially decay to its stationary 
state by i N l = A N . For the SSEP(K N ), A(K N ) = E N {N /2 - 
1 ) = N is the same in every invariant sector of constant par- 
ticle number of the process, except in the one-dimensional 
sectors of n = and n = N, for which there is no gap at all. 
That A(Kpj) = N hints at the fact that the characteristic time 
Ttv scales with the system size as NTn = NA(Kpj)~ l = 1, i.e., 
the interacting particle system relaxes to its stationary state 
after just one step, irrespective of N. This has to do with the 
fact that on Ajy any vertex can be reached from any other one 
through a single jump. We also see that the energy levels are 
monotone decreasing in S, E^(S) > E^(S') whenever S < S', a 
property that has been dubbed "ferromagnetic ordering of en- 
ergy levels" for ferromagnetic SU(2)-invariant quantum spin 
models [16]. There are both examples and counterexamples 
of models holding this property [17]. These results are akin to 
the so-called Aldous' spectral gap conjecture (ca. 1990), ac- 
cording to which the gap of the single particle random walk 
(n = 1) should be equal to the gap of the interchange pro- 
cess (n = AO on any finite simple graph [18]. This conjecture 
spawned some original results in probability and mathemati- 
cal physics, mostly on the last decade, and was proved in the 
most general form only recently through a melange of group- 
theoretical, probabilistic, and combinatorial arguments [19]. 



IV. THE SSEP ON THE COMPLETE BIPARTITE GRAPH 

The complete bipartite graph Kn { ^ is the simple undi- 
rected graph with partitioned vertex set V = Vi U V2 with 
\V{ I = Ni, i = 1 , 2, and V\ n V2 = <2> such that every vertex in V\ 
is connected to every vertex in V2 by a unique edge; see Fig. 1. 



TABLE I. Eigenspectrum of the SSEP(G) on the complete graph A"g. 
The degeneracies g;v(5) are given as (25+ 1) x / 2 (iV,S). Notice 
that Y,sSn(S) = 2 8 , as it should be. 



5 1 2 3 4 
E N {S) 20 18 14 8 
g N {S) 1x 14 3x 28 5x 20 7x79x1 



TABLE II. Characterization of the invariant subspaces Q.„ of the 
SSEP(AT 8 ). The multiplicities of the \S,M = -N/2 + n) states within 
each fl„ are given in the last column as (5 d i/ 2 ( JV ' S )). 





-M dimn n 







1 


(4 1 ) 


1 


8 


(3 7 )(4 1 ) 


2 


28 


(2 20 )(3 7 )(4') 


3 


56 


(1 28 )(2 20 )(3 7 )(4 1 ) 


4 


70 


(0 14 )(1 28 )(2 20 )(3 7 )(4>) 


5 


56 


(1 28 )(2 2 «)(3 7 )(4') 


6 


28 


(2 20 )(3 7 )(4 1 ) 


7 


8 


(3 7 )(4>) 



8 1^ (4 1 ) 



For this graph, 

H = \ E L ( l ~°h ■ ofe) = ^1^2-251 -S 2 , (14) 
i 1 ev 1 ! 2 eV2 L 

where the operators S\ and Si are given by 

Si = \ £ a h , 5 2 = ~ £ G h . (15) 

iieVi Z ( 2 gv 2 

It is a matter of simple algebra to demonstrate that the Si, i — 
1,2, obey Si x Si = iSi, being thus legitimate spin operators. 
The magnitude of the spin Si is 5, = Ni/2. 

Notice that in representing the occupation state of V, by 
a state of 5, indexed by the value of its S-' component 
through the relation n,- = A^,72 +#»,-, m/ = —Ni/2,— Ni/2 + 
1, . . . , +N/2, we have promoted a reduction of the dimen- 
sion of the configuration space associated with V; from 2 ' 
to Ni + l. This reduction comes from lumping equivalent con- 
figurations obtained by permutations of the particles among 
the vertices of Vi into a single representative state. The result 
is that the 2 A '-dimensional original problem can be treated as 
a (A^i + 1)(A^2 + 1) -dimensional problem as far as the determi- 
nation of the eigenspectrum is concerned. If the eigenstates of 
(14) become needed, e.g., to calculate correlation functions or 
block entropies, one must reconstruct them from the original 
2 N states by appropriate combinations of permutations. 

In terms of the total spin S = Si+ Si, we have 2S\ ■ Si = 
S 2 — Sf — S\- In the basis diagonal in the complete set of com- 
muting operators S\, S^S 2 , and S z = S\+Sj, the eigenvalues 
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of (14) are given by 

E NltN2 (Si,S2,S,M) = 

1 (16) 
= + Si (Si + 1) + S 2 (S 2 + 1) - S(S + 1), 

or, in more compact form, by 

E N (S) = (N/2-S)(N/2 + S+l), (17) 

with |Si-S 2 | sC 5 ^ Si + 5 2 and \M\ < 5. The number of 
particles in the system is given by n = N/2+M, as before. 

For each of the min{25i + 1,252 + 1} values of 5, En(S) 
is 25 + 1 degenerate due to its independence on M. Over- 
all, M is in the range — |5i — S 2 \ ^ M ^ Si +52, and for 
any given M we have max{|M|,|5i — 52|} ^5 5 ^ 5i + 52- 
The lowest eigenvalue of (14) lies in the sector of maxi- 
mum 5 = 5i + 52 = N\/2 +N 2 /2 = N/2 — as expected for 
a "ferromagnetic" model — , with a N +1 degeneracy (M = 
-N/2, -N/2+l,..., +N/2) associated with the N + 1 sta- 
tionary states of the process, one within each invariant sector 
of constant number of particles (« = N/2 +M = 0,1,.., ,N). 
The steady states are given by the same |<J>^ («)) as in Eq. (13). 

The spectral gap of the process is given by A^a^ a^) = 
E N {N/2-\)=N + 2, and like the gap of the SSEP&v) is 
the same in every invariant sector of constant particle number. 
It is also clear from Eq. (16) or Eq. (17) that the energy lev- 
els observe the ferromagnetic ordering property mentioned in 
Sec. Ill, namely, E^{S) > En(S') whenever 5 < 5', providing 
yet another example of such systems [16, 17]. 



V. THE SSEP ON MULTIPARTITE GRAPHS 

The cases analyzed so far lead naturally to the SSEP on gen- 
eralized multipartite graphs. In particular, two types of multi- 
partite graphs are of interest: complete mutipartite graphs and 
concatenated (chained) bipartite graphs. Although the SSEP 
on this second type of graphs cannot be solved by elementary 
techniques — actually, some of them cannot be exactly solved 
at all — , they give rise to infinitesimal generators that may ap- 
peal in other modeling circumstances. 



A. The complete multipartite graph 

The complete multipartite graph Kq 1 ... q n is the simple 
undirected graph with partitioned vertex set V — V\ U ■ ■ ■ U Vn 
with \Vi\ =Qi,i=l,... ,N, and V, HV/ = for / + 1 j such that 
every two vertices from different sets Vj and Vj are adjacent. 
When Qi = ■ ■ ■ = Qn = Q, we have the Q-regular complete 
multipartite graph Kg ; see Fig. 2. 

Following our previous approach, we associate to each dis- 
joint subset Vk a spin-6/t/2 operator Sk acting on its own sub- 
space of dimension Qk + l given by 




FIG. 2. The 3-regular complete multipartite graph Kj 3 3 = if| occu- 
pied by n = 5 particles (black circles). Clearly, K\ 1 = A*^ = K^, 

the complete graph considered in section III. 

In terms of these spin operators, the infinitesimal generator of 
the SSEP on the complete multipartite graph Kq { .g w reads 

1=1 (=1 /=i 1=1 

On the Q-regular complete multipartite graph Kg, all 5, are 
equivalent spin-g/2 operators. In this case, taking into ac- 
count that Sf = 56(76+1)' th e infinitesimal generator (19) 
oftheSSEP(^) becomes 

X= \nQ{\nQ+\)- (£a) 2 - (20) 

This operator is formally identical with the Hamiltonian of a 
quantum spin-6/2 Curie-Weiss model and can be analyzed 
along the same lines as the spin-j operator (12) in Sec. III. If 
we define a total spin operator S = Si, then the eigenvalues 
of (20) can be read off immediately as 

E N Q {S,M)= l -NQ( l -NQ+\)-S{S+\), (21) 

with 5 = 5 min ,5 min + \,...,NQ/2 and M = -5,-5 + 
l,...,+5, where 5 mm = 1/2 if Q and are both odd and 
Srnin = otherwise. Within each sector of fixed number of 
particles n — NQ/2 +M, the values of 5 range in the inter- 
val \M\ ^ 5 ^ NQ/2. The degeneracies associated with the 
eigenvalues (21) are given by gw(5) = (25+ 1) x cIq/ 2 (N,S), 
where now the outer multiplicities cIq/ 2 (N,S) determining the 
degeneracies in the 5 values are given by (cf. appendix) 

d Q/2 (N, 5) = b Qj2 (N, S)-b Q/2 (N,S+l), (22) 

where the coefficients bQ/ 2 (N,M) are given by 

b Q/2 (N,M) = 

= Y(-i) k ( N ) (ihQ+^N+M-iQ+m-n 

h WV \QN + M-{Q+\)k )' 

(23) 



Sk= ~ L d k- 



(18) 



where the summation runs over k as long as the summing 
terms are non-null. Equations (22)-(23) are analogous to 
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TABLE III. Characterization of the invariant subspaces £l n for the 
SSEP^j). The dimensionality of Q.„ is given by the coefficient 
b Q / 2 {N,M = \-NQ/2 + n\) given in Eq. (23). The multiplicities of 
the \S,M = —NQ/2 + n) states within each Sl„ are given in the last 
column as (S'WC" 5 )). The additional multiplicities coming from 
the permutation equivalent states within each set Vq., i = l,...,N, 
are not accounted for in this table. 



n= \NQ + M dim fl„ {S d ai l{N ' S) ) 



o 

\J 


1 




1 


4 


(5 3 )(6') 


2 


10 


(4 6 )(5 3 )(6 1 ) 


3 


20 


(3 I0 )(4 6 )(5 3 )(6 1 ) 


4 


31 


(2 11 )(3 10 )(4")(5 3 )(6 1 ) 


5 


40 


(1 9 )(2 ii )(3 io )(4 6 )(5 3 )(6 i ) 


6 


44 


(0 4 )(l 9 )(2 n )(3 10 )(4 6 )(5 3 )(6 


7 


40 


(1 9 )(2 n )(3 io )(4 6 )(5 3 )(6 i ) 


8 


31 


(2 1I )(3 10 )(4 6 )(5 3 )(6 1 ) 


9 


20 


(3 I0 )(4 6 )(5 3 )(6 1 ) 


10 


10 


(4 6 )(5 3 )(6') 


11 


4 


(5 3 )(6') 


12 


1 


(6 1 ) 



Eq. (12) and, indeed, they reduce to it when Q = 1. Co- 
efficient (23) also corresponds to the dimension of the in- 
variant subspace £2m of fixed M = |— NQ/2 + n\; notice that 
dimO^ = dim£2_M- 

The same observations made for the eigenspectrum (10) 
hold here. The zero eigenvalue in Eq. 21) occurs in the sec- 
tor of S = NQ/2 with a NQ + 1-fold degeneracy, correspond- 
ing to the NQ + 1 stationary states of the process, one within 
each subspace of constant particle number. The spectral gap 
A(Kq ) = NQ of the process is also the same within each in- 
variant sector of constant particle number. Finally, it is clear 
from Eq. (21) that the eigenvalues observe the ferromagnetic 
ordering E%(S) > E%(S') if S < S'. In fact, the SSEP(^) 
and the SSEP(X^) differ only by the total spin associated with 
each vertex set, barring the additional degeneracies induced 
by the permutational equivalence of states within each vertex 
set that we briefly discussed in Sec. IV. 




FIG. 3. Concatenated bipartite graphs in a chain. The segment shown 
depicts graphs K42, X/34, and K42 linked together such that every ver- 
tex is linked to every other vertex in the "neighboring" vertices sets. 
In general, the edges set of the composite graph under open bound- 
ary condition is given by E opea = (V Ql x Vq 2 ) U (Vq 2 x Vg 3 ) U • ■ ■ U 
(^Gw x ^G«+i)' un der closed boundary condition the edges set be- 
comes Closed = £open U (V Qf/+1 X V Qi ). 



where each spin operator Si is given as in Eq. (18). Opera- 
tor (24) is equivalent to the Hamiltonian of a one-dimensional 
ferromagnetic Heisenberg model of mixed spins S\ = Qi/2, 
. . . , Sn+i = Qn+\/2, with each bipartite graph Kq u q m of the 
chain corresponding to a "unit cell." 

If all Qi are equal, then a simple spin-wave analysis shows 
that the low-lying eigenspectrum just above the stationary 
state has the form En « 2 sin 2 (ft /./V), with an asymptotic be- 
havior E?t ~ 2n 2 N~ 2 for N > 1 [14]. In this case, the rela- 
tion between the relaxation time scale and the spectral gap be- 
comes T N 2 , typical of diffusive behavior. It is well known 
that the SSEP displays this type of dispersion relation, where 
the dependence on the number of particles (and on the spin 
Q/2) affects only prefactors, not the dependence on A^ 2 [1, 2]. 

If some or all Qi are different, then we have a full-fledged 
mixed-spins operator. It has been proved that the eigenspec- 
trum of the mixed-spins chain (24) displays the ferromag- 
netic ordering of energy levels; actually, this property was first 
demonstrated for quantum spins chains exactly like (24) under 
open boundary conditions [16]. While mixed-spins Hamilto- 
nians of interest in the theory of magnetism are usually anti- 
ferromagnetic, operator (24) is always ferromagnetic [28]. We 
thus expect that a modified spin-wave analysis already suc- 
cesful in the more complicated cases of antiferromagnetic or 
competing interactions [30] shall work even better in the anal- 
ysis of (24). This provides an interesting avenue for further 
investigations. 



B. Concatenated bipartite graphs in a chain 

If we concatenate bipartite graphs, we obtain a graph like 
the one depicted partly in Fig. 3. We shall denote this graph 
as Kq 1 q 2 x Kq 2 q 7> x • • • x Kq n q n+1 . Under periodic bound- 
ary conditions there is an additional link between Kq n ,q n+1 
and Kq^Qt, and in this case we must have Qn+i = Qi, other- 
wise it is an open chain. For this graph, under open boundary 
conditions 5i reads 

^=i£afitfi-2£A-Afi, (24) 
L 1=1 1=1 



VI. SUMMARY AND OUTLOOK 

The investigation of interacting particle systems on graphs 
is an active field of mathematical reasearch [11-13, 19]. In the 
physics literature, however, one is almost always concerned 
with linear chains, hypercubes, and complete graphs, this last 
ones in the guise of mean-field approximations. We showed 
that other types of graphs are amenable to investigation with 
familiar techniques like quantum angular momentum algebra 
and basic group representation theory. The list of graphs that 
can be explored in this way includes star graphs (the star graph 
Sn is just the complete bipartite graph K\ wheel graphs, 
and finite regular trees; see Fig. 4. 
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FIG. 4. Left: The star graph S5. Star graphs Sn are equivalent to 
complete bipartite graphs Kin. Right: The wheel graph W5. 



We avoided employing SU(A0 representations and Young 
tableaux on purpose to keep the exposition elementary. In- 
deed, the study of permutation invariant operators like (8) and 
(20) is more natural by means of the permutation group; see, 
e.g., [20] for background and [21-23] for closely related ap- 
plications. Otherwise, it is clear that interacting spin waves, 
variational states (including matrix product ansatze), and clus- 
ter approximations, among other approaches, could also be 
applied in the investigation of exclusion processes on graphs, 
e.g., to estimate the gap of the mixed-spins operator (24). 

While several works have investigated the quantum Curie- 
Weiss model and related XY and Ising models in transverse 
fields in the context of quantum phase transitions and quan- 
tum information and computation [24-27], very little, if any, 
has been done on the Heisenberg model on complete bipar- 
tite and multipartite graphs. These models deserve a closer 
look, because they are very simple and because complete bi- 
partite and multipartite graphs appear in the discussion of sev- 
eral combinatorial problems that may become relevant in the 
field of quantum information and computation. 
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Appendix: Outer multiplicities in the Clebsch-Gordan series for 
arbitrary spin-s representations of the SU(2) 

The degeneracies of the eigenvalues of the operators (8) and 
(20) partly come from the existence of many possible different 
combinations of the elementary spins summing up to a definite 
value of S. This degeneracy is encoded in the outer multiplic- 
ity d s (N,S) of the irreducible representation (D^ appearing 
in the A^-fold tensor product 



sN 



(A.l) 



Un- 



resolved into the direct sum via repeated application of the 
Clebsch-Gordan series [20] 

2)W ® = 2>(l'-<l) © ©(I^M) © ... © 2>('+Q , (A.2) 



where 5 m ; n = 1/2 if i is half-integer and N is odd and S m m = 
otherwise, N ^ 2 is the number of vertices of the graph, and s 
is the magnitude of the spins involved. 

To calculate d s (N,S), we first notice that in a system of N 
spins s, the number b s (N,M) of states of total magnetization 
M is given by the coefficient of z M in the expansion of (z~ s + 

z~ s+l H \-zf) ■ This recipe stems from the solution of the 

simple combinatorial problem of distributing M things among 
N boxes each supporting a minimum of —s and a maximum of 
+s things [31]. Next we notice that d s (N,S) and the numbers 
b s (N 7 M) are related by 



d s (N,S)=b s (N,S)-b s (N,S+l) 



(A3) 



since in a given subspace of fixed M we have S ^ \M\, so that 
b s (N,S) — b s (N,S + 1) counts just those states with exactly 
total spin S. Finally, an explicit expression for b s (N, S) can be 
obtained from its defining expansion, 



(z- s + z- s+l + --- + z 



s\N . 



+sN 
M=-sN 



= z (1+zH \-z ) = z 



2sN 

~ sN 52 c 2s+i{N,M)z M , 

M=0 

(A.4) 



such that b s (N,M) = c^+i (N,sN + M). It is clear from 
Eq. (A.4) that b s {N,M) = b s {N,-M). The coefficients 
c 2s+i(v) are known as generalized (or extended) binomial 
coefficients of order 2s + 1, and reduce to the standard bino- 



mial coeficients when s = \ [31, 32], 



b l/2 (N,M)=c 2 (N,lN+M) 



N 
N + M 



(A.5) 



compare Eqs. (A. 3) and (A.5) with Eq. (12). If we put z = 1 
in Eq. (A.4) we obtain Y,M b s(N,M) = (2s + l) N , as required. 

It turns out that generalized binomial coefficients can be 
written in terms of standard binomial coeficients [31, 32]. The 
resultant expression for b s (N,M) is 



N\ f(s + l)N+M-(2s+l)k-l 
sN + M-(2s+l)lc 



b s (N,M) = Y,(- 

(A.6) 

where the summation runs over k as long as the summing 
terms are non-null. Both the upper and the lower terms in 
the second binomial coefficient above are integer, even if Af 
is odd and s is half-integer, because then M will necessar- 
ily be half-integer. Equations (A. 3), (A.4), and (A.6) to- 
gether solve the Clebsch-Gordan series decomposition (A.l) 
for arbitrary spin-s representations of the SU(2). For exam- 
ple, we find that the irreducible representation ©( 29 / 2 ) appears 
^3/2(15,29/2) = 167895 times in the Clebsch-Gordan series 
for [2)( 3 / 2 )]® 15 after computing just 10 terms in Eq. (A.6) 
for each coefficient ^3/2(15, 29/2) and £3/2(15, 31/2) of 
Eq. (A. 3), a modest effort for an otherwise laborious calcu- 
lation. In fact, it is possible to combine (by means of ele- 
mentary binomial identities) the two coefficients in Eq. (A. 3) 
into a single factor, further reducing the calculational effort. 
Of course, in practice d s (N,S) is calculated on a computer, 
except perhaps in small cases (say, for s ^ 1 and Af ^ 10). 
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